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The data on arrival directions of ultra-high energy cosmic rays (UHECRs) detected with the Yakutsk array are 
analyzed. The work is induced by the recent claim of the Pierre Auger collaboration for the significant correlation 
found between UHECRs and positions of nearby Active Galactic Nuclei (AGN) on the celestial sphere; and no 
correlation the HiRes collaboration stands for. Conflicting data of four giant arrays concern possible extragalactic 
sources of UHECRs and appeal to the profound analysis and to the future data from the Telescope Array/Northern 
Auger Observatory. 



1. Introduction 

Detection of extensive air showers (EASs) of 
cosmic rays and the evaluation of energies, masses 
and arrival directions of primary particles is a 
crucial technique to search for UHECR sources. 
The localization of the sources is complicated due 
to charged particle trajectories distorted in un- 
known (extra)galactic magnetic fields. Only at 
the highest energies a deflection of protons, the 
most probable particles of cosmic rays, is less than 
or comparable to the angular resolution of the gi- 
ant arrays. 

The Pierre Auger Observatory (PAO) collabo- 
ration has recently analyzed a sample consisting 
of 81 EASs with energies above 40 EcV (= 4x 10 19 
eV) detected from January 1, 2004 to August 31, 
2007 1 . The authors used a part of the data 
(to May 27, 2006) in order to determine the pa- 
rameters resulting in the maximum correlation of 
UHECR arrival directions with AGN. Then, the 
second part of the data was used to confirm the 
hypothesis obtained. 

As a result, the observed UHECR arrival di- 
rections are found to be anisotropic, there is a 
significant correlation of EASs at energies above 
56 EeV within an angle of i/j = 3.1° with AGN 
from catalog [2] located at distances z < O.Ol&Q 



from the Earth. In the second part of the data 
(from May 27, 2006), 8 of 13 EASs correlate with 
AGN under the same conditions that have been 
found for the first part, while the number of ex- 
pected coincidences is 2.7 in the isotropic case, 
with the chance probability P = 1.7 x 10 -3 . 

This result is confirmed by the Yakutsk array 
data [3], while the HiRes data demonstrate no 
significant correlation with AGN 0]. 

Gorbunov et al. stated contrary to the PAO 
hypothesis: the conclusion that the bulk of UHE- 
CRs are protons originating in nearby AGN can 
be rejected at 99% CL. Instead, they attribute 
PAO observational data to the existence of a 
bright source in the direction of the Centaurus Su- 
percluster [5] . Another interpretation of the data 
was proposed by Wibig and Wolfendalc, namely, 
that cosmic rays are nuclei with In A = 2.2. ± 0.8 
generated in radio galaxies [6 . 

In this article the situation is re-visited using 
available data of the giant surface arrayt@ and 
taking into account the HiRes result. 



x 75 Mpc, assuming H = 71kms 1 Mpc 



2 e.g. PAO is a surface array where the energy is estimated 
via the fluorescence light; at the Yakutsk array the energy 
estimation is based on the Cherenkov light measurements. 
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2. The Yakutsk array and experime 
data used in analysis 

The Yakutsk array geographical coordiu 
are 61. 7° V, 129A°E, 100 m above sea leveH 
present it consists of 58 ground-based and 6 
derground scintillation detector stations to i 
sure charged particles (electrons and muons) 
48 detectors of the air Cherenkov light. Du 
more than 30 years of lifetime the array has 1 
re-configured several times, the total area cov 
by detectors was maximal about 1990 (S <- 
km 2 ), now it is S ~ 10 km 2 . During the w 
observation period approximately 10 6 showers of 
the primary energy above 10 15 eV have been de- 
tected. 

In this work a sample of the data published by 
Pravdin et al. [7] (51 EASs with energies above 
40 EeV and zenith angles below 60°) is analyzed. 
The angular resolution error is less than 5° for 
these showers. The energy estimation method is 
based on the total flux measurement of the air 
Cherenkov light and the number of electrons and 
muons at observation level [8 9 J . The energy esti- 
mation error is about 30% and 50% for the show- 
ers with axes inside the array area and in the ef- 
fective region outside, respectively. In some cases 
a data sample was extended down to E t hr = 1 
EeV in order to trace the energy dependence. 

Additionally, the EAS events detected with 
AG AS A (58 EASs, E > 40 EeV, [TU]) and PAO 
(27 EASs, E > 57 EeV, [TT]) are used in analy- 
sis. Arrival directions of UHECRs under consid- 
eration are illustrated in Fig. [T] in supergalactic 
coordinates using Hammer- Aitoff projection, to- 
gether with celestial positions of AGN from the 
catalog [2], with redshifts z < 0.018. 

3. Exposure of sky zones to the surface ar- 
rays 

Due to the array acceptance area depending 
on zenith angle and energy because of the shower 
attenuation in the atmosphere as wel as the geo- 
metric factor, the calculation of a celestial region 
exposure to the surface array in the diurnal cycle 
is a non-trivial task. 
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Figure 1. UHECR arrival directions in su- 
pergalactic coordinates. Observational data: 
AG AS A (circles); PAO (squares); the Yakutsk 
array (triangles); AGN from Veron's catalog 
(crosses). Galactic plane is indicated by dots. 



A simple method to calculate the exposure is 
the Monte Carlo algorithm implementation [12] , 
In this reversal approach one can take N random 
points distributed isotropicalljQ in the horizontal 
system, transform to equatorial angles and than 
count the points hit within the given celestial re- 
gion. In the limit N — > oo the exact exposure of 
the region is evaluated. 

In Fig.[5]the resultant relative exposures of the 
right ascension zones (a £ (0°, 360°), 3°j < 6 < 
3°(j + 1), j = —30, ..,30) to three surface arrays 
are shown; N = 10 7 . Fields of view are com- 
plementary in the case of the Yakutsk array and 
PAO, while AG AS A is observing a region in low 
latitudes where is a depression in the consolidated 
fields of other two arrays. 

4. Searching for correlation with the large 
scale structures 

The distribution of luminous matter is highly 
anisotropic in the universe, notably, there are 
two planar structures - Milky Way disk within 
the Galaxy and the so-called supergalactic plane 
(SGP) [13] in the galaxy distribution. There may 



3 1020 g/cm 2 . 



4 f(9, <f>) oc sin 29 at energies above 10 EeV. 
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Figure 2. Sky coverage by the surface 
arrays: AGASA (35.8°iV, 6 max = 45°); 
PAO (32.2°S, B max = 60°) and Yakutsk 



(61.7° AT, ( 



60°). E > 10 EeV. 



Figure 3. Galactic plane enhancement differ- 
ence between observed, i?, and expected in the 
isotropic case, Rq, values, d = 10°. Vertical bars 
are statistical errors given by Eq. [2l horizontal 
bars indicate energy bins. 



be a correlation between these structures and 
UHECR arrival directions, if the concentration 
of potential sources (e.g. SNRs, AGN) is propor- 
tional to the density of matter. 

In this section, the Yakutsk array data are re- 
analyzed in comparison with the data from other 
arrays in order to verify this possibility. We have 
used the (super)galactic plane enhancement pa- 
rameter, R, which is given by 

n(|6| < d)-n(\b\ > d) 



R 



n(\b\ < d)+n(\b\ > d)' 



(1) 



where b is the (super)galactic latitude; d is the 
plane border. A statistical error of the parameter 
is 



SR 



' n(\b\ < d)n(\b\ > d) 
(n(\b\ < d)+n(\b\ > d)) 1 



(2) 



A galactic plane enhancement (GPE) parame- 
ter is estimated using the available data of three 
giant arrays: AGASA, PAO and Yakutsk. In the 
latter case EAS events are selected with energies 
above 1 EeV to reveal the energy dependence of 
the excess flux, if any. The results are shown in 
Fig. [3] An expected parameter is calculated by 
the Monte Carlo method (N = 10 6 ) using the 
algorithm described in the previous section. 

There is no excess flux from the galactic disk 
in the whole energy range examined. The local 



excess in the interval (10, 20) EeV is insignifi- 
cant. Angular dimension of the possible sources 
disk is uncertain, so we have scanned the range 
d 6 (0°,30°). None of d values resulted in the 
significant GPE. 

Another possibility for the anisotropy in 
UHECR arrival directions is connected with the 
SGP, as was hypothesized by Stanev et al. [T4] . 
In order to reveal the possible excess we have cal- 
culated the enhancement parameter, SGPE, and 
its statistical error repeating the procedure above 
with definitions (|1I2|) for a supergalactic latitude. 

A scan over the SGP angular dimension d £ 
(0°,30°) revealed a maximum enhancement ef- 
fect, SGPE, at d = 10°. The difference with 
regard to the 'isotropic' value expected for the 
arrays is shown in Fig. [5] as a function of energy. 

Again, there is no statistically significant en- 
hancement from SGP, in all the data in any en- 
ergy bin. But it seems to be a systematic increase 
of the SGPE in the Yakutsk array data which is 
supported by the data from AGASA and PAO 
above 40 EeV: all the observational enhancement 
parameters are greater than those expected in the 
isotropic case. The possible utility of this result 
consists in expected evidence of the excess flux 
at the highest energies from future arrays, or in- 
creased number of EAS events from the PAO ar- 
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Table 1 

The number of coincidences, Nhiti m arrival di- 
rections of N EASs above it^r with AGN. N iso 
and P are an expected number and a chance prob- 
ability for the isotropic distribution. 
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24 
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Figure 4. Supcrgalactic plane enhancement pa- 
rameter, d = 10°. 



ray. 

5. Searching for possible extragalactic 
sources of UHECRs 

In this section we are going to verify the PAO 
claim with independent datasets provided by 
AGASA and the Yakutsk arrays. In contrast with 
the previous analysis [5] , we have used as the null 
hypothesis the maximum correlation parameters 
given by the whole PAO dataset scanned above 40 
EeV, rather than predictor-confirmer scheme di- 
viding a sample into two parts. These parameters 
are: a correlation within angle ip — 3.2° between 
AGN at z < 0.017 and 20 UHECRs out of the 
total 27 with energies above E t hr = 57 EeV [TT] . 

Observed and expected-for-isotropy numbers of 
UHECRs correlated with AGN are calculated for 
the Yakutsk and AGASA EAS events applying 
the sky exposure (Fig. \2§ to arrays. Resultant 
numbers and chance probabilities are given in Ta- 
ble Q] together with the PAO result. 

The null hypothesis is rejected in AGASA case, 
while the Yakutsk array data confirm the result of 
the PAO collaboration, but at a lower significance 
level. An analysis of the HiRes data has resulted 
in no significant correlation with AGN [4] , so the 
pro and con arguments are in equality. Appar- 
ently, we have to wait for the future data in this 
case, too. 

Due to the different fields of view, as well as the 



energy /arrival angles estimation procedure, the 
'optimal' correlation parameters can be different 
for the arrays. For this reason, a scan in the en- 
ergy (E > 40 EeV), redshift (0.001 < z < 0.03) 
and angular distance (1° < ijj < 6°) of the 
Yakutsk array data is performed to determine the 
maximum ratio of the difference in the observed 
number of coincidences and the number expected 
for the isotropic case to the standard deviation. 
The maximum ratio appears to be reached for 22 
EASs with energies above 60 EeV, 12 of which 
arrive within tp = 3° of the AGN (while an ex- 
pected number is 4.1) at the distance from the 
Earth less than z = The chance proba- 

bility is P = 2 x 10~ 4 . However, it is necessary 
to apply a penalty factor to the probability due 
to a posteriori selection of the parameters in this 
case. 

Figure [5] shows the ratio of the observed num- 
ber of coincidences in UHECR arrival directions 
(within ip — 3°) with quasars, Lacertae, and AGN 
from the catalog [5] to the number of random 
coincidences expected for the isotropic distribu- 
tion. The active nuclei are chosen at the distances 
z < 0.015, the quasars are taken with redshifts 
z < 0.3, and BL Lacs are selected with the lumi- 
nosities m < 18, as in [15] . It is found that there 
is no significant deviation of the coincidences from 
isotropic expectation in the case of Lacertae and 
quasars. A variation in the redshift for quasars 
does not reveal any significant deviation. No sig- 
nificant correlation is also found for HP objects 
and BL+HP objects from the catalog [2J. The 
significant correlation with AGN is shown in Fig- 
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Figure 5. Ratio of the coincidences of UHECR 
arrival directions with extragalactic objects to the 
number of random coincidences expected for the 
isotropic distribution versus the threshold energy 
of the particles, E > Eth r \ the Yakutsk array 
data. The statistical error bars are also shown. 



ure as was found scanning parameters above. 

It is interesting to know whether all classes 
of the objects belonging to AGN are sources or, 
for instance, only Seifert galaxies generate UHE- 
CRs as was suggested by Uryson [16] . To answer 
this question, one can analyze a correlation of 
UHECR arrival directions with AGN classes sep- 
arately. Figure [5] shows the results. AGN objects 
are divided into four classes along the proposal of 
VeronCetty and Veron [2j: i) SI, Seifert galax- 
ies of the first type with broad Balmer lines; ii) 
S2, Seifert galaxies of the second type; iii) S3, so- 
called LINERs (low-ionization nuclear emission- 
line regions), which are galaxies with weak nu- 
clear emission lines; and iv) H2, galaxies whose 
spectrum of nuclear emission lines is similar to 
that of nebulae ionized by hot stars. The redshift 
boundaries giving the maximum correlation with 
the Yakutsk array data are selected for each class 
of the objects: z < 0.015 for SI and S3, z < 0.016 
for S2, and z < 0.024 for H2. 

Correlations of Seifert galaxies of the first 
type and LINERs with UHECRs do not exceed 
'isotropic' expectation. Only S2 and H2 objects 
correlate with UHECRs (with the maxima at 



Figure 6. Ratio of UHECR correlation with dif- 
ferent classes of AGN to the expected number. 



E > 50 and E > 70 EeV, respectively). The ex- 
cess in the observed number of coincidences over 
the number expected for the isotropic distribu- 
tion is 4.7<t and 3.7<t in the maxima for S2 and 
H2. Therefore, possible sources of UHECRs are 
Seifert galaxies of the second type and/or H2 ob- 
jects at the distances less than 100 Mpc. 

6. Correlations in the redshift intervals 

Another part of the PAO hypothesis is that 
UHECRs correlate with AGN at distances z < 
0.018 due to the Greisen-Zatsepin-Kuzmin ef- 
fect [17] , which strongly suppress the flux of cos- 
mic rays with energies E > 60 EeV from cos- 
mological distances. In order to verify this ef- 
fect with the Yakutsk array data, the ratio of 
the observed number of coincidences of cosmic 
rays (E > 60 EeV) with AGN in various red- 
shift intervals to the number of random coinci- 
dences expected for the isotropic distribution is 
used (Fig. [?])■ Indeed, a significant correlation of 
UHECRs with AGN is found only in the near- 
est bin z G (0.001,0.015). In all other redshift 
bins, the observed number of coincidences is equal 
(within errors) to the number expected in the 
isotropic case. Hence, this can be considered as 
one of the independent evidences of the Greisen- 
Zatsepin-Kuzmin effect. 
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Table 2 

Correction factors to energy scales of the EAS 
array pairs, Re, averaged in the region E > 1 
EeV. 
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Figure 7. Ratio of the coincidences of UHECR 
arrival directions with AGN to the expected num- 
ber in various redshift bins, z. The vertical bars 
are statistical errors. The boundaries of z bins 
(shown by the horizontal bars) are 0.001, 0.015, 
0.025, 0.05, 0.075, 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.7, 
0.9, 1.5, and 5.4. 



7. Corrections to the measured UHECR 
intensity and energy 

Correlations with AGN are found in the PAO 
and Yakutsk array data almost within the same 
angular spots (3.2° and 3°, respectively). At the 
same time, an approximate coincidence in the 
threshold energies at which the maximum corre- 
lation is observed by the two arrays is acciden- 
tal. A comparison of the energy spectra measured 
with different arrays infer the existence of system- 
atic differences between cosmic ray energies esti- 
mated. Additionally, another correction factor to 
the measured intensity of UHECRs is caused by 
the instrumental errors and fluctuations of the 
shower parameter^]. The latter factor, which is 
derived to be [15] : 

<t 2 k 2 

Rj = exp{—^-), (3) 

where a is RMS deviation; n is the integral en- 
ergy spectrum index, should be applied before 
any comparison of the energy spectra observed 
at different EAS arrays. 

6 for a given primary particle energy 



While the true energy correction factors are un- 
known for arrays, we can cross calibrate the en- 
ergy estimation methods adjusting arbitrary fac- 
tors, Re, for the pairs of spectra measured, which 
converge them properly. The resulting spread 
of factors elucidates a confidence interval for the 
UHECR energy estimated. 

Table [2] demonstrates the variety of correction 
factors to the energy estimation methods used. 

To illustrate the result of corrections Rj x Re 
applied to estimated energies and intensities, the 
measured spectra are shown in Fig. [8j Energy 
scale factors are used here from the second col- 
umn of Table [2l although any other column may 
be used as well. 

The UHECR propagation model results are 
presented by the two examples: Bahcall and Wax- 
man model (shown by dots) [15], and Berezinsky 
et al's (solid curve) [20] . 

In this context, the actual energy thresholds 
of the AGN correlation effect in the case of the 
Yakutsk array and PAO data seem badly different 
(Table [3D. 



Table 3 

Energy thresholds where the maximum correla- 
tion of UHECRs and AGN is found. 
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